30th International Cosmic Ray Conference 




ICRC? 



Effect of muon-nuclear inelastic scattering on high-energy atmospheric muon 
spectrum at large depth underwater 

S. I. SlNEGOVSKY 1 , A. MlSAKI 2 , K. S. LOKHTIN 1 AND N. TAKAHASHI 3 . 
^Irkutsk State University, Irkutsk, Russia 

2 Advanced Research Institute for Science and Technology, Waseda University, Tokyo, Japan 

3 Department of Advanced Physics, Hirosaki University, Hirosaki, Japan 
misakiakeo @air. ocn, ne.jp 

Abstract: The energy spectra of hadron cascade showers produced by the cosmic ray muons travelling 
through water as well as the muon energy spectra underwater at the depth up to 4 km are calculated with 
two models of muon inelastic scattering on nuclei, the recent hybrid model (two-component, 2C) and 
the well-known generalized vector-meson-dominance model for the comparison. The 2C model involves 
photonuclear interactions at low and moderate virtualities as well as the hard scattering including the 
weak neutral current processes. For the muon scattering off nuclei substantial nuclear effects, shadowing, 
nuclear binding and Fermi motion of nucleons are taken into account. It is shown that deep underwater 
muon energy spectrum calculated with the 2C model are noticeably distorted at energies above 100 TeV 
as compared to that obtained with the GVMD model. 



Introduction 

The muon inelastic scattering on nuclei contributes 
noticeably to the total energy loss of cosmic rays 
muons. The influence of this interaction on the 
shape of ultra-high energy muon spectra at the 
great depth of a rock/water is still unknown in 
detail. Of interest is also to estimate the num- 
ber of cascade showers produced by very high- 
energy muons through inelastic interactions with 
nuclei, and to study the influence of this pro- 
cess on the energy spectra of cosmic -ray muons 
in water at the depths of the underwater/ice neu- 
trino telescopes - NT200+, AMANDA, IceCube, 
ANTARES, NESTOR, NEMO and KM3NeT. 

In this work, we calculate energy spectra of hadron 
cascade showers produced by cosmic -ray (atmo- 
spheric) muons in water due to inelastic scattering 
on nuclei, as well as the integral energy spectra of 
atmospheric muons in water at depths up to 4 km. 
Calculations are performed with two models: the 
hybrid model of inelastic scattering of leptons on 
nuclei [1,2] and, for a comparison, the known gen- 
eralized vector-meson-dominance (GVMD) model 



of photonuclear muon interactions by Bezrukov 
and Bugaev [3]. 

Muon-nucleus inelastic scattering 

The hybrid two-component (2C) model [1,2] for 
inelastic scattering of high-energy charged leptons 
on nuclei involves photonuclear interactions at low 
and moderate Q 2 as well as the deep inelastic scat- 
tering processes at high Q 2 values. For the vir- 
tually < Q 2 < 5 GeV 2 the Regge based 
parametrization [4] for the electromagnetic struc- 
ture function is applied and the muon-nucleon 
inelastic scattering cross section at Q 2 < 5 GeV 2 
is computed with the formula 



da 
dQ 2 dy 



2(1 + R) 




2m 2 



(1) 



where the ratio R = a^ /ax is taken into account 
according to Ref. [5], y = v/E is the fraction of 
muon energy transfered to the hadron system, Q 2 
and x — Q 2 / (2MEy) are the Bjorken variables. 
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Table 1 : Comparison of calculations of the muon energy loss due to inelastic scattering in standard rock. 



E, 




bn(E), 


10" b cm 2 - g- L 






GeV 


[1,2] 


[14] 


[13] 


[15] 


[16] 




(2C7) 










10 5 


0.62 


0.60 


0.68 


0.70 


0.70 


10 6 


0.82 


0.80 


0.88 


1.08 


1.00 


10 8 


1.53 


1.60 




2.25 


2.50 


10 9 


2.16 


2.18 




3.10 


4.00 



In the range Q 2 > 5 GeV 2 the cross section of 
/ziV-scattering can be written in the form 
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where we put R = Q 2 jv 2 , that is equivalent to 
the Callan-Gross relation, F 2 = 2xF 1 . Signs "±" 
stand for fi ± . In Eq. (2) used notations are: 

F?° = - g^zFf + (g» 2 + gf)v 2 z F 2 z , 

F 3 NC = -9^zFf + 2g^zF 3 Z ; (3) 
G F M 2 Q 2 



2V2ira M 2 + Q 2 ' 

-i + 2sin 2 ^, = ~. 



(4) 
(5) 



The structure functions (SFs) F 2 , F% represent 
the weak neutral current (NC) contribution, F^ Z , 
F^ z are taking into account the electromagnetic 
and weak current interference. The nucleon SFs, 
F 2 7 , F^ z , F 2 Z , F^ 2 , F 3 z , are defined in the quark- 
parton picture through the parton distributions (see 
[6, 7]). For the range of Q 2 > 6 GeV 2 the 
electroweak nucleon SFs are computed with the 
CTEQ6 [8] and MRST [9] sets of the parton dis- 
tributions. Linear fits for the nucleon SFs are used 
in the range 5 < Q 2 < 6 GeV 2 . 

Nuclear modifications of the nucleon SFs due to 
coherent and incoherent effects [10, 11] are taken 
into account according to Ref. [12] (see also [13]) 
by the factor r A (x, Q 2 ) = F^/F^ that is the ratio 
of the structure function per nucleon of a nucleus 
with mass number A and the isospin averaged nu- 



cleon structure function F 2 JV . Thus, 
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The contribution of muon-nuclear interactions to 
the continous energy loss is given by integral 
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where Nq — Na/A is the number of nuclei per 
gramme of matter. 

Table 1 presents the muon energy loss b n due to 
muon-nucleus interactions in rock, computed for 
various models: values of b n obtained for the 2C 
model are given in the second column, the rest are 
results predicted in Ref s. [13, 14, 15, 16]. One may 
see that the muon energy loss b n (E) at E > 10 6 
GeV obtained in Refs. [1,2, 14] differ apparently 
from the predictions [15, 16]. 

The energy dependence of the inelastic scattering 
energy loss for muons traveling through water may 
be parametrized in the energy range 10 2 — 10 9 GeV 

as 



b n (E) = c + cirj + c 2 7? 2 + c 3 ri 3 + c^rf 



(8) 



where 77 = log 10 (£'/l GeV) and coeffecients Cj 
(in units of 10 -6 cm 2 g _1 ) are 

c = 1.06416, ci = -0.64629, c 2 = 0.20394, 
c 3 = -0.02465, c 4 = 0.00130. (9) 

Figure 1 shows the energy loss b n (E) for the in- 
elastic scattering of muons in water calculated with 
the 2C model (solid curve) and the GVMD one 
(dashed). 
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Figure 1: The muon energy loss due to muon- 
nucleus interactions in water. 

Muon induced hadron showers in water 

The number of hadron showers with energies 
above ui per cm 2 per second per steradian, gen- 
erated in water column Ah = h 2 — hi through 
muon-nucleus interactions at energy above E, may 
be defined as 

S n (u,E,Ah,6) =N J dh J de D^(e,h,8)x 



dy (10) 



dy 



Here D M {E,h, 



the 
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muon flux 
9 is zenith 

angle, y max = 1 — rn^/E. Ratios of the shower 
spectra computed with the 2C model to those ob- 
tained with the GVMD one are shown in figure 2. 
At the muon energy E = 10 TeV numbers of the 
muon-induced hadron showers, calculated with 
two models, the 2C and GVMD, differ slightly 
(about 10%). However the discrepancy between 
the models grows with increasing energy: for 
E = 100 TeV, it is as great as ~ 30%, and for 
E = 10 5 TeV, the result obtained for the 2C model 
in the range of the large energy transfer exceeds 
by a factor of ~ 2.5 the corresponding result for 
the model of photonuclear interaction [3]. For 
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Figure 2: The ratio of the energy spectra of hadron 
showers computed with the 2C model to those ob- 
tained with use of the GVMD one. 

energies E > 10 3 TeV, the number of showers cal- 
culated with the 2C model exceeds that obtained 
with use of the GVMD model by 20-50% even 
in the region of small energy loss (y ~ 0.1). For 
the catastrophic energy loss (y > 0.5), the number 
of showers, obtained with use of the 2C model, 
exceeds that of the GVMD prediction by factor 
about 2. 



Atmospheric muon spectra underwater 

The solution of the muon transport equation [17] 
allows to obtain the differential energy spectra 
D^{E,h : 9) at large depth in homogeneous me- 
dia and therefore to compute the flux of cosmic ray 
muons with energy above specified one: 



N^{E,h,t 



deD fJi (e, h,( 



(11) 



Figure 3 shows the ratio 

N^ c {E,h)/N^ VMD (E,h) of the near verti- 
cal muon flux underwater computed at depth 1-4 
km for two models of muon-nuclear scattering, 
the 2C model [2] and GVMD one [3]. The effect 
is rather noticeable at the E > 10 TeV, especially 
for the depth 3-4 km. For h = 4 km this ratio 
decreases to 0.75 at E — 10 3 TeV. Thus a sizeable 
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Figure 3: The ratio of muon fluxes underwater cal- 
culated with the 2C model to those obtained with 
the GVMD one. 

increase of the muon inelastic scattering cross 
section may result in an appreciable decrease 
of the deep underwater muon flux as compared 
to that obtained [18, 19] for the GVMD model. 
It should be noted that this result refers only to 
the atmospheric conventional (n, K) muons. As 
concerns muons produced in charmed particle de- 
cays (prompt muons), which become presumably 
dominant at E > 10 5 GeV (see e. g. [18, 19]), the 
role of the muon-nucleus inelastic scattering needs 
further study. 

Conclusions 

Evidently the increase of the cross section of in- 
elastic muon scattering in matter, while leading to 
diminished cosmic -ray muon flux deep underwa- 
ter, results in growing efficiency of muon registra- 
tion. This last factor is positive for neutrino astron- 
omy since neutrino-induced muons may yield the 
signal from astrophysical high-energy muon neu- 
trinos. 
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